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ABSTRACT 

In the unification scheme, radio quasars and FR II radio galaxies come from the same parent population, but 
viewed at different angles. Based on the Comptonization models for the 7-ray emission from active galactic 
nuclei (AGNs), we estimate the number of radio quasars and FR II radio galaxies to be detected by the Gamma- 
Ray Large Area Space Telescope (GLAST) using the luminosity function (LF) of their parent population derived 
from the flat-spectrum radio quasar (FSRQ) LF. We find that ~ 1200 radio quasars will be detected by GLAST, 
if the soft seed photons for Comptonization come from the regions outside the jets. We also consider the 
synchrotron self-Comptonization (SSC) model, and find it unlikely to be responsible for 7-ray emission from 
radio quasars. We find that no FR II radio galaxies will be detected by GLAST. Our results show that most radio 
AGNs to be detected by GLAST will be FSRQs (~99% for the external Comptonization model, EC model), 
while the remainder (~1%) will be steep-spectrum radio quasars (SSRQs). This implies that FSRQs will still 
be good candidates for identifying 7-ray AGNs even for the GLAST sources. The contribution of all radio 
quasars and FR II radio galaxies to the extragalactic 7-ray background (EGRB) is calculated, which accounts 
for -30% of the EGRB. 

Subject headings: galaxies: active — galaxies: jets — accretion, accretion disks — radio continuum: galaxies 



1. INTRODUCTION 

The third catalog of 7-ray AGNs detected by the Ener- 
getic Gamma-Ray Experiment Telescope (EGRET) on the 
Compton Gamma-Ray Observatory (CGRO) includes ~ 80 
high-confidence identific ations of blazars (e.g., Hartm arTet al.l 
1999; Matt ox et al.ll2001l) . GLAST has higher sensitivity than 
EGRET, and much more blazars are expected to be detected 
after its launch. Many workers have predict ed the statistic 
properties of b l azars in the G LAST era (e.g.. ISteckerl 1 1 999t 
lPadovanill2007t iDerm er 2007). One method is to extrapolate 
the observed 7-ray luminosity distribution of EG RET blazars 
to ob tain a 7-ray luminosity function (LF) (Chiang et al. 
1995). An alternative method is to assume some correlation 
between 7-ray emission and the emission in other bands to 
model the undetected 7-ray blazars, in which the larger sam- 
ples in other bands provide useful clues to such researches 
(e.g..lPadovani et alJ ll993t IStecker et alJl!993b lDermeiil2007t 
lPadovanill2007l) . The previous works on the EGRB showed 
that about ~ 25% to ~ 100% of th e EGRB can be at- 
tributed to the unresolved blazars (e.g., Padovani et al. |ll993t 
Chiang et al]l!995t IStecker & Salamon ll 19961: iMuc ke & Pohll 
2000; Narumoto & Totanill2006l) ~ 

Comptonization is widely believed to be responsible for the 
7-ray emission from the blazars detected by EGRET, which 
can be classified into two categories: the EC models and 
SSC model, a ccording to the origin of the soft seed pho- 
tons (see e.g., iBottcherl 120071 for a review and references 
therein). The space density and evolution of the parent pop- 
ulation of blazars, together with the Lorentz factor distribu- 
tion of the jets, are crucial for understanding the properties 
of 7-ray emitting blazars. In almost all previous works, the 
models of blazars are rather simplified. In this Letter, we 
derive the parent radio LF of radio quasars/FR lis from the 
FSRQ LF to investigate the statistic properties of 7-ray emit- 
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ting quasars to be detected by GLAST. The cosmological pa- 
rameters f^M = 0.3, Ha = 0.7, and Ho = 70 km s" 1 Mpc" 1 have 
been adopted in this Letter. 

2. EC AND SSC MODELS FOR 7-RAY AGNS 

In the EC models, the observed 7-ray emission from the 
relativistic jet is closely related to its observed radio emission 
(see Eq. 26 in lDermer etaD 19971) . 
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where u* is the soft seed photon energy density measured in 
the stationary source frame, mh is the magnetic energy density 
in the jet, /i b s = cos# b s (6* b s is the direction of the jet motion 
with respect to the line of sight), Sj = [7j(l -/3j^ bs)] _1 for a jet 
moving at /3jC, and vU^ is the observed radio luminosity of 
the jet. The energy distribution of the nonthermal electrons in 
the jet is assumed to be n e oc j^ P - 

In the EC models, the soft photons may originate from the 
accretion d isks, the broad-line regions ( BLRs), or/and the dust 
tori (e.g., [ Ghisellini & Mada ul 119961 ; iGeorgano poulos et alj 
I200U IDermer & Schlickeiserll2002l) . It was argued that the 
contribution from the accretion disks is not important, be- 
cause the 7-ray emitting region is far away from the disk 
and in the jet comoving frame, the energy density of photons 
from the disk is deboosted by the relativistic je t moving away 
from the black hole (e.g., ISikora etail Il994t IDermer et alj 
11997b . It is well known that the BLR size Rblr oc L^j LR 
where a BhR ~ 0.5 -0.7 ( e.g., iKaspi et a l. 2000; Ben tz et alj 
120061) . iBentz et all ([2006) found that a B LR is 0.518 subtract- 
ing the contribution from the host galaxy starlight to Lboi, 
which is consistent with q;blr — 0.5 expected from the photo- 
ionization model if all BLRs have similar physical properties. 
The inner radius of the dust torus is roughly at th e dust evap- 
oration radius: R mn oc L°^ dNetzer & Laorlfl993l) . The pho- 
ton energy density u* oc L/R 2 , where L = Lblr and R = Rblr 
for the BLR photons, and L = Lir and R ~ R inn for the dust 
torus case. The irradiated infrared luminosity of the dust torus 



2 



Lir oc Lboi, if the opening angle o f the toru s does not vary 
much for individual sources fe.p.. lCac3l2005h . Thus, the en- 
ergy density of the soft photons from the BLRs/dust tori is 
roughly universal for most sources. We rewrite Eq. ([T]i as 



where i= 1, re, and 
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where the normalization Cec is related with u* /mh (see Eq. 
[T]). For t he SSC model, the o bserved 7-ray luminosity is (see 
Eq. 28 in lDermer et aII[T997h 
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where Cssc is the normalization. 



3. THE PARENT LF OF FSRQS 

In the unification scheme, FRSQs, SSRQs, and FR II galax- 
ies come from the same parent population, but viewed at dif- 
ferent angles. Some previous authors have derived the par- 
ent LF using differ ent approaches (see ICara & Listeri 120071 : 
iLiu & Zhand 120071) . however, they have not compared the 
n umber density of b l azars with that of radio galaxies. 

iPadovani & Urrvl d 19921) derived the radio LFs of FSRQs 
and FR II galaxies from a sample of radio-loud AGNs. They 
considered a two-component model, in which the total lumi- 
nosity L„ t is the sum of an unbeamed part £„ u and a jet 
luminosity 

(4) 
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They used a variety of observational features to constrain 
the ratio /(= C v JC VyU ). They found that a constant / ~ 
4.5 x 10~ 3 can succ essfully explain the observations (see 
IPadovani &Urryll 19921 for the details). Thus, FSRQs should 
satisfy <5j > <5j, m j n ~ 6.45, as their core dominance parameters 
R = L v JC v . a > 1 is required, where an average core spec- 
tral index a ra d = -0.1 is adopted. They further assumed that 
the probability distribution of the Lorentz factors for the jets 
is P(7j) = C7j G , between 7^1 = 5 and 7^2 = 40. Their de- 
rived FSRQ LF is consistent with the beaming model and 
the LF of FR II galaxies, pro vided G = -2.3 is adopted. Re- 
cently, Pa dovani et al.l d2007l) derived a FSRQ LF based on 
the deep X-ray radio blazar survey (DXRBS) in the same 
way, which extends to l ower luminosity than that derived by 
IPadovani & Urrvl dl992l) . 

The sources in the parent population may be observed as 
FR lis, when thei r jets are oriented at ang les # bs > 40° to the 
line of sight (e.g., Padov ani & Urrvl 19921) . The sources in this 
parent population with o t, s < 40° and 5j < <5j. m ; n will appear 
as SSRQs. The LFs of FSRQs, SSRQs, or FR II galaxies 
can be reproduced with this parent radio LF, if the probability 
distribution of the Lorentz factors P(7j) is supplied. The LF of 
FSRQs 0FSRQ(£iAj) can be derived from the LF of the parent 
population 4>(C u j) by 

7j,2 1 

f f dj £v i 

0FSRQ(i 1 /j,z)= / f(7j) d 7j / (KCj^-j^-d^obs, (5) 

where the orientations of the jets of the parent population 
are assumed to be isotropically distributed, and only those 
with ^ > M ™ n (7j) = (7Amm- l)/(7f- l) I/2 <5j.mm are FS- 
RQs, which is required by FSRQs having R > 1. Equation 
([5]) can be re-written as 
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The coefficient e« can be calculated with Eq. (f7]) by using Eq. 
© and adopting the integral limits /i bs(7j) to satisfy £*j — 
ACvj/2 < C u j < C'l ■ + ACvj/2. Solving a set of n linear 
algebraic equations (|6| numerically, the parent LF 4>(C v ^z) 
ca n be calculated from th e LF of FSRQs ^fsrq^i/jiZ) given 
bv lPadovani et all d2007l) . 



4. NUMBER OF GLAST QUASARS 

Using this derived parent radio LF, we can calculate the ob- 
served 7-ray LF 4> 1 for FSRQs, SSRQs and FR II galaxies 
based on either EC or SSC models: 
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where dC v j/dL^ SSC can be derived from Eqs. (f2]l-(|4li. This 
7-ray LF is not limited to blazars, as their parent radio LF 
is adopted in the calculations. Adopting the conditions for 
different sources (i.e., <5j > <5j, m in for FSRQs; 6* b s > 40° for FR 
II galaxies; and the remainder are SSRQs), the numbers of the 
FSRQs/SSRQs/FR II galaxies with f v>1 > ff* as functions of 
redshift z can be calculated by 



cW EC / ssc (z) 
dz 



= / ^f /SSC (Liy sSC ,z)dLiy ssc , (9) 



either for EC or SSC models, respectively. 

Sixty three 7-ray emitters were identified as blazars at 
high c onfidence with measured redshifts in Hartm arTet alJ 
(1999). More 7-ray sources were identified as blazars af- 
terwa rds dMattoxet alJl200U ISowards-Emmerd et al.l l2003l 
l2004t ISguera et al.ll2004h~ We collect all these 7-ray blazars 
identified at high confidence, which leads to 64 quasars and 
16 BL Lac objects with measured redshifts. We note that the 
measured fluxes with photon energy greater than 100 MeV 
> 5 x 10~ 8 photons cm" 2 for all EGRET blazars. This can be 
translated to vf^ GREI ~8x 10" 12 erg s _1 cm " 2 at 100 MeV 
assum ing a mean photon spectral i ndex of 2 dHartman et alJ 
1999), which corresponds to p = 3 dDermer et al.lll997h . 

We calculate the number counts of 7-ray quasars with 7- 
ray flux densities greater than /™ n at 100 MeV using Eq. 
(0, based on the different 7-ray radiative models (EC or 
SSC models). We tune the values of the parameters Cec or 
Cssc to let the total number of FSRQs derived with Eq. (0 
equal to that of the FSRQs detected by EGRET. We find that 
Cec = 4.42 x 10~ 3 or Cssc = 10.62 are required to reproduce 
64 FSRQs detected by EGRET for the EC or SSC models, 
respectively. Based on the derived values of Cec or Cssc, the 
number counts of 7-ray emitting quasars/FR lis to be detected 
by GLAST can be predicted by adopting different flux density 
limits /„. m i n with Eq. (O (see Fig. [T]and Table 1). 
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FIG. 1. — Upper panel: the redshift distributions of the EGRET blazars 
(the black solid line: quasars; the dashed line: all EGRET blazars; the dotted 
line: BL Lac objects). The color lines represent our model calculations for 
64 quasars. The red line represents the model prediction based on the EC 
model, while the blue line represents the result for SSC model. Lower panel: 
the calculated redshift distributions of 7-ray quasars detected with different 
flux density limits. The red lines represent the calculations based on the EC 
model, while the blue lines are for the SSC model. The dotted lines repre- 
sent our model calculations for EGRET quasars. The dashed lines represent 
the redshift distribution of the 7-ray quasars with > 0. 1/,^™ E1 at 100 MeV, 
while the solid lines are for the sources with > 1 /30/J?^ ET . 

The total numbers of 7-ray emitting radio quasars/FR lis as 
functions of sensitivity are plotted in Fig. [2] About 1200 7- 
ray radio quasars will be detected by GLAST based on the EC 
model, if its s ensitivity is 30 times higher than that of EGRET 
at 100 MeV dGehrels & Michelsonlll999h . Our calculations 
show that no FR II radio galaxies will be detected by GLAST 
as 7-ray emitters either for EC or SSC models. We find that 
almost all 7-ray quasars (~99%) to be detected by GLAST 
will be FSRQs for the EC model, and the remainder (~1%) 
will be SSRQs. For the SSC model, - 1800 quasars will be 
detected by GLAST, of which -80% will be FSRQs (see Fig. 
[2]). We use the derived 7-ray LF (Eq. [8]) to calculate the con- 
tribution of all radio quasars/FR lis to the EGRB (listed in 
Table 1). 




limit limit 

FIG. 2. — The total number of the 7-ray quasars with flux densities greater 
than /ii m it at 100 MeV. The red lines are calculated for the EC model, while 
the blue lines are for the SSC model. The dashed lines represent the 7-ray 
FSRQs, while the colored dotted lines represent the 7-ray SSRQs. The dash- 
dotted line represents the FR II radio galaxies. The black dotted line repre- 
sents the GLAST sensitivity at 100 MeV. 



5. DISCUSSION 

We find that the redshifts of almost all EGRET BL Lac ob- 
jects are < 1, which implies that BL Lac objects may have dif- 
ferent space density and evolutionary behaviors from quasars. 
The LF of BL Lac ob jects was derived from the DXRBS by 
iPadovani et alj d2007l) . however, the results for BL Lac ob- 
jects are more uncertain than those for FSRQs, because of the 
small number statistics and —30% of them having no redshift. 
In this work, we use a parent radio LF of radio quasars/FR 
II galaxies derived from the FSRQ LF. The derived redshift 
distributions of 7-ray quasars are similar for different models 
(EC or SSC), which are roughly consistent with that of the 
EGRET quasars. 

It was suggested that the 7-ray radiative mechanisms are 
different for quasars and BL Lac objects, i.e., the EC mech- 
anism may be respon sible for quasars, while th e SSC is for 
BL Lac objects (e.g.. iDondi & Ghisellinilll995h . If the EC 
mechanism is indeed responsible for 7-ray quasars, the pre- 
dicted 7-ray quasars to be detected by GLAST will be — 1200. 
Our results are roughly consistent with the estimate given by 
iDermeri (|2007) based on a simplified blazar model. The SSC 
model predicts a simple relation between 7-ray luminosity 
and radio luminosity of the jets (see Eq. 0. The EGRET 
flux limit ^/"egret at MeV can be converted to a radio 
flux density limit /„, ra d ~ 10 Jy at 5GHz by using Eq. (f3]). 
This means that all EGRET quasars should have their radio 
flux densities higher than — 10 Jy, which is obviously incon- 
sistent with most EGRET quasars h aving f VfS & > 1 Jy (e.g., 
IStecker etalll 19931: IZhou et all [l 9971) . Thus, the SSC model 
is unlikely to be responsible for EGRET quasars, unless the 
physical properties of the jets are significantly different for 
individual sources, i.e., the values of Cssc f° r most sources 
deviate significantly from a constant value. 

Our results show that no FR II galaxies will be detected by 
GLAST. Most GLAST quasars will be FSRQs (-99% for the 
EC model), which implies that FSRQs will still be good can- 
didates for identifying the 7-ray sources even for the GLAST 
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TABLE 1 

The summary of different models 



Cec 


Cssc N(>^™)a 


N(> 1/10*/™) 


N(> 1/30*/™) 


-'EGRB 


4.42 xl(T 3 


64 

10.62 64 


491 

696 


1203 
1806 


0.32 
0.30 



11 The number of FSRQs with > vffgjg ^ 8 X 10" 
above 100 MeV in units of 10~ 5 photons cm -2 s _1 sr - 

sources. At l east two FR I g alaxies have been detected by 
EGRET (e.g.,[Padovani 2007), and more FR I galaxies were 
predi cted to be detected by GLAST (e.g., iGhisellini et alj 
2005). In the unification scheme, FR I galaxies are BL Lac 
objects with misaligned jets. The 7-ray FR I galaxies/BL 
Lac objects are beyond the scope of this work. The EGRET 
quasars identified at high confidence (64 sources) have ra- 
dio flux densities > 1 Jy and are about 20% of the total FS- 
R Qs (~ 300) above the same flux density limit (see Fig. 6 



Padovani et al.l [20071) . Assuming the radio flux density 
limit of the GLAST quasars to be ~ 30 times lower than 
the EGRET limit, i.e., ~ 30 mJy, the total all-sky number of 
the FSRQs above this limit is about 20000. So, about 4000 
FSRQs will be detected by GLAST, if the same percentage 
(~ 20%) is adopted as the EGRET blazars. This rough esti- 
mate is about a factor of two higher than our model calcula- 
tions. Considering that some unidentified 7-ray sources are 
likely to be blazars, our estimates on 7-ray quasars to be de- 
tected by GLAST are only lower limits. Both the duty cycle 
in the 7-ray band and identification rate may affect the detec- 
tion rate, and we implicitly assume them to be similar to those 



s cm " at 100 MeV. for the photons with energy 

of EGRET blazars. 

The EGRB integrated above 100 MeV was determined 
to be 1.45(± 0.05)x 10~ 5 photons c m" 2 s" 1 sr" 1 from the 
EGRET data (ISreekumar et al.l [19981) . IStrong et all (120041) 
used a new model of the Galactic background, and obtained a 
slightly smaller value of the EGRB, 1.14(±0.12) x 10" 5 pho- 
tons cm" 2 s" 1 sr" 1 . We sum up the 7-ray emission from all 
radio quasars/FR lis with our derived 7-ray LF and find that 
they contribute ^30% of the EGRB (see Table 1), which are 
only lower limits, because our calculations are limited to ra- 
dio quasars/FR II galaxies. The BL Lac objects/FR Is must 
contribute some part to the EGRB. The detailed calculation 
of their contribution to the EGRB is beyond the scope of this 
Letter. 
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